Markers of fertility in reproductive microbiomes of male and female endangered black-footed ferrets (Mustela nigripes).

Reproductive microbiomes contribute to reproductive health and success in humans. Yet data on reproductive microbiomes, and links to fertility, are absent for most animal species. Characterizing these links is pertinent to endangered species, such as black-footed ferrets (Mustela nigripes), whose populations show reproductive dysfunction and rely on ex-situ conservation husbandry. To understand microbial contributions to animal reproductive success, we used 16S rRNA amplicon sequencing to characterize male (prepuce) and female (vaginal) microbiomes of 59 black-footed ferrets at two ex-situ facilities and in the wild. We analyzed variation in microbiome structure according to markers of fertility such as numbers of viable and non-viable offspring (females) and sperm concentration (males). Ferret vaginal microbiomes showed lower inter-individual variation compared to prepuce microbiomes. In both sexes, wild ferrets harbored potential soil bacteria, perhaps reflecting their fossorial behavior and exposure to natural soil microbiomes. Vaginal microbiomes of ex-situ females that produced non-viable litters had greater phylogenetic diversity and distinct composition compared to other females. In males, sperm concentration correlated with varying abundances of bacterial taxa (e.g., Lactobacillus), mirroring results in humans and highlighting intriguing dynamics. Characterizing reproductive microbiomes across host species is foundational for understanding microbial biomarkers of reproductive success and for augmenting conservation husbandry.

The gut microbiomes of Channel Island foxes and island spotted skunks exhibit fine-scale differentiation across host species and island populations.

California's Channel Islands are home to two endemic mammalian carnivores: island foxes (Urocyon littoralis) and island spotted skunks (Spilogale gracilis amphiala). Although it is rare for two insular terrestrial carnivores to coexist, these known competitors persist on both Santa Cruz Island and Santa Rosa Island. We hypothesized that examination of their gut microbial communities would provide insight into the factors that enable this coexistence, as microbial symbionts often reflect host evolutionary history and contemporary ecology. Using rectal swabs collected from island foxes and island spotted skunks sampled across both islands, we generated 16S rRNA amplicon sequencing data to characterize their gut microbiomes. While island foxes and island spotted skunks both harbored the core mammalian microbiome, host species explained the largest proportion of variation in the dataset. We further identified intraspecific variation between island populations, with greater differentiation observed between more specialist island spotted skunk populations compared to more generalist island fox populations. This pattern may reflect differences in resource utilization following fine-scale niche differentiation. It may further reflect evolutionary differences regarding the timing of intraspecific separation. Considered together, this study contributes to the growing catalog of wildlife microbiome studies, with important implications for understanding how eco-evolutionary processes enable the coexistence of terrestrial carnivores-and their microbiomes-in island environments.

Epigenetic changes to gene pathways linked to male fertility in ex situ black-footed ferrets.

Environmental variation can influence the reproductive success of species managed under human care and in the wild, yet the mechanisms underlying this phenomenon remain largely mysterious. Molecular mechanisms such as epigenetic modifiers are important in mediating the timing and progression of reproduction in humans and model organisms, but few studies have linked epigenetic variation to reproductive fitness in wildlife. Here, we investigated epigenetic variation in black-footed ferrets (Mustela nigripes), an endangered North American mammal reliant on ex situ management for survival and persistence in the wild. Despite similar levels of genetic diversity in human-managed and wild-born populations, individuals in ex situ facilities exhibit reproductive problems, such as poor sperm quality. Differences across these settings suggest that an environmentally driven decline in reproductive capacity may be occurring in this species. We examined the role of DNA methylation, one well-studied epigenetic modifier, in this emergent condition. We leveraged blood, testes, and semen samples from male black-footed ferrets bred in ex situ facilities and found tissue-type specificity in DNA methylation across the genome, although 1360 Gene Ontology terms associated with male average litter size shared functions across tissues. We then constructed gene networks of differentially methylated genomic sites associated with three different reproductive phenotypes to explore the putative biological impact of variation in DNA methylation. Sperm gene networks associated with average litter size and sperm count were functionally enriched for candidate genes involved in reproduction, development, and its regulation through transcriptional repression. We propose that DNA methylation plays an important role in regulating these reproductive phenotypes, thereby impacting the fertility of male ex situ individuals. Our results provide information into how DNA methylation may function in the alteration of reproductive pathways and phenotypes in artificial environments. These findings provide early insights to conservation hurdles faced in the protection of this rare species.

Are rabid raccoons (Procyon lotor) ready for the rapture? Determining the geographic origin of rabies virus-infected raccoons using RADcapture and microhaplotypes.

North America is recognized for the exceptional richness of rabies virus (RV) wildlife reservoir species. Management of RV is accomplished through vaccination targeting mesocarnivore reservoir populations, such as the raccoon (Procyon lotor) in Eastern North America. Raccoons are a common generalist species, and populations may reach high densities in developed areas, which can result in contact with humans and pets with potential exposures to the raccoon variant of RV throughout the eastern United States. Understanding the spatial movement of RV by raccoon populations is important for monitoring and refining strategies supporting the landscape-level control and local elimination of this lethal zoonosis. We developed a high-throughput genotyping panel for raccoons based on hundreds of microhaplotypes to identify population structure and genetic diversity relevant to rabies management programs. Throughout the eastern United States, we identified hierarchical population genetic structure with clusters that were connected through isolation-by-distance. We also illustrate that this genotyping approach can be used to support real-time management priorities by identifying the geographic origin of a rabid raccoon that was collected in an area of the United States that had been raccoon RV-free for 8 years. The results from this study and the utility of the microhaplotype panel and genotyping method will provide managers with information on raccoon ecology that can be incorporated into future management decisions.

Gut bacterial composition shows sex-specific shifts during breeding season in ex situ managed black-footed ferrets.

The gut microbiome of mammals engages in a dynamic relationship with the body and contributes to numerous physiological processes integral to overall health. Understanding the factors shaping animal-associated bacterial communities is therefore paramount to the maintenance and management in ex situ wildlife populations. Here, we characterized the gut microbiome of 48 endangered black-footed ferrets (Mustela nigripes) housed at Smithsonian's National Zoo and Conservation Biology Institute (Front Royal, Virginia, USA). We collected longitudinal fecal samples from males and females across two distinct reproductive seasons to consider the role of host sex and reproductive physiology in shaping bacterial communities, as measured using 16S rRNA amplicon sequencing. Within each sex, gut microbial composition differed between breeding and non-breeding seasons, with five bacterial taxa emerging as differentially abundant. Between sexes, female and male microbiomes were similar during non-breeding season but significantly different during breeding season, which may result from sex-specific physiological changes associated with breeding. Finally, we found low overall diversity consistent with other mammalian carnivores alongside high relative abundances of potentially pathogenic microbes such as Clostridium, Escherichia, Paeniclostridium, and (to a lesser degree) Enterococcus - all of which have been associated with gastrointestinal or reproductive distress in mammalian hosts, including black-footed ferrets. We recommend further study of these microbes and possible therapeutic interventions to promote more balanced microbial communities. These results have important implications for ex situ management practices that can improve the gut microbial health and long-term viability of black-footed ferrets.

Novel genomic resources contribute to the systematics of threatened arboreal deer mice of the genus Habromys Hooper & Musser, 1964 (Cricetidae, Neotominae) within a neotomine-peromyscine phylogeny.

The Crested-tailed deer mouse, Habromyslophurus, is one of seven arboreal species within the genus Habromys. Species of this genus are monotypic, relatively rare, and occur in low densities. Their geographical distribution is highly fragmented due to being restricted to montane cloud forest in Mesoamerica and they are of conservation concern. All Habromys species are endemic to Mexico, except H.lophurus, which is also distributed in Guatemala and El Salvador. In this study, we obtained and characterized the first mitogenome and several thousand nuclear ultraconserved elements (UCEs) of H.lophurus to determine its phylogenetic position within neotomine-peromyscine mice. Its mitogenome sequence (16,509 bp) is only the second complete mitogenome obtained for this poorly known genus. We also obtained the first nuclear genomic data for H.lophurus, including 3,654 UCE loci, as well as a partial mitogenome of H.simulatus (6,349 bp), and 2,186 UCE for the outgroup Holochilussciureus. Phylogenetic analyses that included our newly generated genomic data coupled with previously published data from other neotomine-peromyscine mice confirm the placement of H.lophurus, H.simulatus, and H.ixtlani within a highly supported clade. The Habromys clade was nested within a clade that also contains members of the genus Peromyscus and provides further support for the hypothesis of the paraphyly of Peromyscus. These genomic resources will contribute to future phylogenomic studies that aim to further elucidate the evolutionary history of this rare and critically endangered genus of rodents.

Chromosome-length genome assembly and karyotype of the endangered black-footed ferret (Mustela nigripes).

The black-footed ferret (Mustela nigripes) narrowly avoided extinction to become an oft-cited example of the benefits of intensive management, research, and collaboration to save a species through ex situ conservation breeding and reintroduction into its former range. However, the species remains at risk due to possible inbreeding, disease susceptibility, and multiple fertility challenges. Here, we report the de novo genome assembly of a male black-footed ferret generated through a combination of linked-read sequencing, optical mapping, and Hi-C proximity ligation. In addition, we report the karyotype for this species, which was used to anchor and assign chromosome numbers to the chromosome-length scaffolds. The draft assembly was ~2.5 Gb in length, with 95.6% of it anchored to 19 chromosome-length scaffolds, corresponding to the 2n = 38 chromosomes revealed by the karyotype. The assembly has contig and scaffold N50 values of 148.8 kbp and 145.4 Mbp, respectively, and is up to 96% complete based on BUSCO analyses. Annotation of the assembly, including evidence from RNA-seq data, identified 21,406 protein-coding genes and a repeat content of 37.35%. Phylogenomic analyses indicated that the black-footed ferret diverged from the European polecat/domestic ferret lineage 1.6 million yr ago. This assembly will enable research on the conservation genomics of black-footed ferrets and thereby aid in the further restoration of this endangered species.

Genomic resources for Asian (Elephas maximus) and African savannah elephant (Loxodonta africana) conservation and health research.

We provide novel genomic resources to help understand the genomic traits involved in elephant health and to aid conservation efforts. We sequence 11 elephant genomes (5 African savannah, 6 Asian) from North American zoos, including 9 de novo assemblies. We estimate elephant germline mutation rates and reconstruct demographic histories. Finally, we provide an in-solution capture assay to genotype Asian elephants. This assay is suitable for analyzing degraded museum and noninvasive samples, such as feces and hair. The elephant genomic resources we present here should allow for more detailed and uniform studies in the future to aid elephant conservation efforts and disease research.

First description of the mitochondrial genomes of the Central American brocket deer Mazama temama (Kerr, 1792) and the Yucatán Peninsula brocket deer Odocoileus pandora Merriam, 1901.

BACKGROUND: The Central American (Mazama temama) and the Yucatán Peninsula brocket deer (Odocoileus pandora) are deer species with cryptic habits, and little is known about their biology. Odocoileus pandora is listed as Vulnerable on the 2015 IUCN Red List of Threatened Species, while M. temama is considered Data Deficient; however, it currently faces a decreasing population trend. METHODS AND RESULTS: We assembled the complete mitochondrial genome for two M. temama specimens and one complete and one partial for O. pandora from Illumina 150 bp paired-end reads. The mitogenomes of M. temama and O. pandora have a length of 16,479-16,480 and 16,419 bp, respectively, AT-biased; they consist of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and one non-coding control region, most of them follow a transcription direction in the heavy strand of the molecule. The mitochondrial genome of O. pandora shows some particularities compared to other deer species, like a shorter control region of 987-990 bp and a cytochrome b gene with a length of 1,143 bp. Our phylogenetic analyses confirm the close affinity of M. temama to South American M. americana and the nested position of the genus Odocoileus, including O. pandora, into the genus Mazama. CONCLUSIONS: Here, we described for the first time the complete mitochondrial genome for these two species. While our study provides additional information about the taxonomic status of the northern neotropical brocket deer, further research is needed to solve the complicated taxonomy of neotropical deer.

First mitochondrial genome of the Amazonian marsh rat Holochilus sciureus Wagner 1842 (Rodentia, Cricetidae, Sigmodontinae).

The Amazonian marsh rat, Holochilus sciureus, is a member of the subfamily Sigmodontinae, the second-largest subfamily of muroid rodents, with 410 species and ca. 84 genera in 12 tribes. This semiaquatic rodent is distributed in South America and is of great economic and epidemiological importance. In this study, we obtained the first mitochondrial genome of the genus Holochilus obtained from a tissue sample associated with a museum voucher specimen. The generated mitogenome sequence of H. sciureus is 16,358 bp length. It comprises a control region and a conserved set of 37 genes encoding for 2 rRNA genes, 22 tRNA genes and 13 protein-coding genes. We conducted a phylogenetic analysis that included H. sciureus and the only five other published mitochondrial genomes of this poorly studied subfamily of rodents.

Old specimens for old branches: Assessing effects of sample age in resolving a rapid Neotropical radiation of squirrels.

Ultraconserved Elements (UCEs) have been useful to resolve challenging phylogenies of non-model clades, unpuzzling long-conflicted relationships in key branches of the Tree of Life at both deep and shallow levels. UCEs are often reliably recovered from historical samples, unlocking a vast number of preserved natural history specimens for analysis. However, the extent to which sample age and preservation method impact UCE recovery as well as downstream inferences remains unclear. Furthermore, there is an ongoing debate on how to curate, filter, and properly analyze UCE data when locus recovery is uneven across sample age and quality. In the present study we address these questions with an empirical dataset composed of over 3800 UCE loci from 219 historical and modern samples of Sciuridae, a globally distributed and ecologically important family of rodents. We provide a genome-scale phylogeny of two squirrel subfamilies (Sciurillinae and Sciurinae: Sciurini) and investigate their placement within Sciuridae. For historical specimens, recovery of UCE loci and mean length per locus were inversely related to sample age; deeper sequencing improved the number of UCE loci recovered but not locus length. Most of our phylogenetic inferences-performed on six datasets with alternative data-filtering strategies, and using three distinct optimality criteria-resulted in distinct topologies. Datasets containing more loci (40% and 50% taxa representativeness matrices) yielded more concordant topologies and higher support values than strictly filtered datasets (60% matrices) particularly with IQ-Tree and SVDquartets, while filtering based on information content provided better topological resolution for inferences with the coalescent gene-tree based approach in ASTRAL-III. We resolved deep relationships in Sciuridae (including among the five currently recognized subfamilies) and relationships among the deepest branches of Sciurini, but conflicting relationships remain at both genus- and species-levels for the rapid Neotropical tree squirrel radiation. Our results suggest that phylogenomic consensus can be difficult and heavily influenced by the age of available samples and the filtering steps used to optimize dataset properties.

Maintenance of local adaptation despite gene flow in a coastal songbird.

Adaptation to local environments is common in widespread species and the basis of ecological speciation. The song sparrow (Melospiza melodia) is a widespread, polytypic passerine that occurs in shrubland habitats throughout North America. We examined the population structure of two parapatric subspecies that inhabit different environments: the Atlantic song sparrow (M. m. atlantica), a coastal specialist, and the eastern song sparrow (M. m. melodia), a shrubland generalist. These populations lacked clear mitochondrial population structure, yet coastal birds formed a distinct nuclear genetic cluster. We found weak overall genomic differentiation between these subspecies, suggesting either recent divergence, extensive gene flow, or a combination thereof. There was a steep genetic cline at the transition to coastal habitats, consistent with isolation by environment, not isolation by distance. A phenotype under divergent selection, bill size, varied with the amount of coastal ancestry in transitional areas, but larger bill size was maintained in coastal habitats regardless of ancestry, further supporting a role for selection in the maintenance of these subspecies. Demographic modeling suggested a divergence history of limited gene flow followed by secondary contact, which has emerged as a common theme in adaptive divergence across taxa.

Phylogenomic systematics of the spotted skunks (Carnivora, Mephitidae, Spilogale): Additional species diversity and Pleistocene climate change as a major driver of diversification.

Four species of spotted skunks (Carnivora, Mephitidae, Spilogale) are currently recognized: Spilogale angustifrons, S. gracilis, S. putorius, and S. pygmaea. Understanding species boundaries within this group is critical for effective conservation given that regional populations or subspecies (e.g., S. p. interrupta) have experienced significant population declines. Further, there may be currently unrecognized diversity within this genus as some taxa (e.g., S. angustifrons) and geographic regions (e.g., Central America) never have been assessed using DNA sequence data. We analyzed species limits and diversification patterns in spotted skunks using multilocus nuclear (ultraconserved elements) and mitochondrial (whole mitogenomes and single gene analysis) data sets from broad geographic sampling representing all currently recognized species and subspecies. We found a high degree of genetic divergence among Spilogale that reflects seven distinct species and eight unique mitochondrial lineages. Initial divergence between S. pygmaea and all other Spilogale occurred in the Early Pliocene (~ 5.0 million years ago). Subsequent diversification of the remaining Spilogale into an "eastern" and a "western" lineage occurred during the Early Pleistocene (~1.5 million years ago). These two lineages experienced temporally coincident patterns of diversification at ~0.66 and ~0.35 million years ago into two and ultimately three distinct evolutionary units, respectively. Diversification was confined almost entirely within the Pleistocene during a timeframe characterized by alternating glacial-interglacial cycles, with the origin of this diversity occurring in northeastern Mexico and the southwestern United States of America. Mitochondrial-nuclear discordance was recovered across three lineages in geographic regions consistent with secondary contact, including a distinct mitochondrial lineage confined to the Sonoran Desert. Our results have direct consequences for conservation of threatened populations, or species, as well as for our understanding of the evolution of delayed implantation in this enigmatic group of small carnivores.

An efficient method for simultaneous species, individual, and sex identification via in-solution single nucleotide polymorphism capture from low-quality scat samples.

Understanding predator population dynamics is important for conservation management because of the critical roles predators play within ecosystems. Noninvasive genetic sampling methods are useful for the study of predators like canids that can be difficult to capture or directly observe. Here, we introduce the FAECES* method (Fast and Accurate Enrichment of Canid Excrement for Species* and other analyses) which expands the toolbox for canid researchers and conservationists by using in-solution hybridization sequence capture to produce single nucleotide polymorphism (SNP) genotypes for multiple canid species from scat-derived DNA using a single enrichment. We designed a set of hybridization probes to genotype both coyotes (Canis latrans) and kit foxes (Vulpes macrotis) at hundreds of polymorphic SNP loci and we tested the probes on both tissues and field-collected scat samples. We enriched and genotyped by sequencing 52 coyote and 70 kit fox scats collected in and around a conservation easement in the Nevada Mojave Desert. We demonstrate that the FAECES* method produces genotypes capable of differentiating coyotes and kit foxes, identifying individuals and their sex, and estimating genetic diversity and effective population sizes, even using highly degraded, low-quantity DNA extracted from scat. We found that the study area harbours a large and diverse population of kit foxes and a relatively smaller population of coyotes. By replicating our methods in the future, conservationists can assess the impacts of management decisions on canid populations. The method can also be adapted and applied more broadly to enrich and sequence multiple loci from any species of interest using scat or other noninvasive genetic samples.

Mitogenome and phylogenetic analyses support rapid diversification among species groups of small-eared shrews genus Cryptotis (Mammalia: Eulipotyphla: Soricidae).

The small-eared shrew genus Cryptotis is the third largest in the family Soricidae and occurs in North, Central, and northern South America. In Mexico and Central and South America, most species inhabit geographically isolated moist, montane habitats at middle and high elevations in a typical sky-island pattern. The 49 recognized species have been partitioned into as many as six species groups based on morphological and molecular phylogenetic studies. The relationships among these species groups are poorly resolved, and their evolutionary histories, including migration patterns and locomotor adaptations, remain unclear. Herein, we provide a new phylogeny incorporating complete mitochondrial genomes (mitogenomes) and supermatrix approach. We compared different evolutionary scenarios using approximately unbiased (AU), Kishino-Hasegawa (KH), and Shimodaira-Hasegawa (SH) statistical tests. The phylogenetic hypothesis based on mitogenomes revealed novel relationships supporting a basal position for the Cryptotis parvus-group in the genus, and a close relationship between C. gracilis and one clade of the C. thomasi-group. The former relationship is consistent with the least derived humerus morphology and northern distribution of the species. The latter relationship implies multiple migrations between Central and South America. The lack of fine resolution for the species group relationships may be due partly to the lack of taxon sampling. In contrast, multi-approach analyses suggest that the unresolved relationships may be a result of rapid diversification during the early stages of Cryptotis evolution.

First mitochondrial genome of the Caucasian squirrel Sciurus anomalus (Rodentia, Sciuridae).

The Caucasian Squirrel, Sciurus anomalus, is the only representative of the Sciuridae family in the Eastern Mediterranean region. In this study, the mitochondrial genome of the Sciurus anomalus species was generated, and we investigate its phylogenetic position within the Sciuridae family. The generated mitogenome sequence is 16,234 bp. It is composed of a control region and a conserved set of 37 genes containing 13 protein-coding genes, 22 tRNA genes and 2 rRNA genes.

First mitochondrial genome of the marbled polecat Vormela peregusna (Carnivora, Mustelidae).

The marbled polecat, Vormela peregusna, is one of the least studied species in the Mustelidae family, especially with regard to phylogeography and genetic diversity. In this study, we determined the mitochondrial genome sequence of V. peregusna and investigated its position within the Mustelidae phylogeny. The generated mitogenome is 15,982 bp in length; it consists of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region.

First mitochondrial genome of the Egyptian mongoose Herpestes ichneumon (Carnivora, Herpestidae).

The Egyptian mongoose, Herpestes ichneumon, is the only extant mongoose in Europe, with populations still distributed in Africa and the Middle East. In this study, we present the first mitochondrial genome sequence of Herpestes ichneumon and we investigate its phylogenetic position within Feliformia suborder. The resultant mitogenome sequence is 16,775 bps, composed of a conserved set of 37 genes containing 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a control region. Our results represent a valuable resource for further phylogeographical studies.

Identification of novel bacterial biomarkers to detect bird scavenging by invasive rats.

Rapid advances in genomic tools for use in ecological contexts and non-model systems allow unprecedented insight into interactions that occur beyond direct observation. We developed an approach that couples microbial forensics with molecular dietary analysis to identify species interactions and scavenging by invasive rats on native and introduced birds in Hawaii. First, we characterized bacterial signatures of bird carcass decay by conducting 16S rRNA high-throughput sequencing on chicken (Gallus gallus domesticus) tissues collected over an 11-day decomposition study in natural Hawaiian habitats. Second, we determined if field-collected invasive black rats (Rattus rattus; n = 51, stomach and fecal samples) had consumed birds using molecular diet analysis with two independent PCR assays (mitochondrial Cytochrome Oxidase I and Cytochrome b genes) and Sanger sequencing. Third, we characterized the gut microbiome of the same rats using 16S rRNA high-throughput sequencing and identified 15 bacterial taxa that were (a) detected only in rats that consumed birds (n = 20/51) and (b) were indicative of decaying tissue in the chicken decomposition experiment. We found that 18% of rats (n = 9/51) likely consumed birds as carrion by the presence of bacterial biomarkers of decayed tissue in their gut microbiome. One species of native bird (Myadestes obscurus) and three introduced bird species (Lophura leucomelanos, Meleagris gallopavo, Zosterops japonicus) were detected in the rats' diets, with individuals from these species (except L. nycthemera) likely consumed through scavenging. Bacterial biomarkers of bird carcass decay can persist through rat digestion and may serve as biomarkers of scavenging. Our approach can be used to reveal trophic interactions that are challenging to measure through direct observation.

First DNA sequence reference library for mammals and plants of the Eastern Mediterranean Region.

The Mediterranean region is identified as one of the world's 36 biodiversity hotspots, with the Earth's most biologically rich yet threatened areas. Lebanon is a hub for Eastern Mediterranean Region (EMR) biodiversity with 9116 characterized plant and animal species (4486 fauna and 4630 flora). Using DNA barcoding as a tool has become crucial in the accurate identification of species in multiple contexts. It can also complement species morphological descriptions, which will add to our understanding of the biodiversity and richness of ecosystems and benefit conservation projects for endangered and endemic species. In this study, we create the first reference library of standard DNA markers for mammals and plants in the EMR, with a focus on endemic and endangered species. Plant leaves were collected from different nature reserves in Mount Lebanon, and mammal samples were obtained from taxidermized museum specimens or road kills. We generated the 12S rRNA sequences of 18 mammal species from 6 orders and 13 different families. We also obtained the trnL and rbcL barcode sequences of 52 plant species from 24 different families. Twenty-five plant species and two mammal species included in this study were sequenced for the first time using these markers.